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Two problems are discussed here. The first one is the 0.4 dex discrepancy between the ^Li abun- 
dance derived from the spectra of metal-poor halo stars on the one hand, and from Big Bang 
nucleosynthesis, based on the cosmological parameters constrained by the Wilkinson Microwave 
Anisotropy Probe (WMAP) measurements, on the other hand. Lithium, indeed, can be depleted 
in the convection zone of unevolved stars, by a combination of diffusion and slow mixing with 
the hotter layers below the convection zone, where ^Li is destroyed by the ^Li(p,a)^He reaction. 
The understanding of the hydrodynamics of this crucial zone near the bottom of the convective 
envelope in dwarfs or turn-off stars of solar metallicity has recently made enormous progress with 
the inclusion of internal gravity waves. However, similar work for metal-poor stars is still lack- 
ing. So in spite of several investigations claiming to have partly explained the above mentioned 
^Li discrepancy, it is not yet clear whether or not the depletion occurring in the metal-poor stars 
themselves is adequate to produce a ^Li plateau. The second problem concerns the difficulties 
met in accounting for the large amount of ^Li recently found in metal-poor halo stars (Asplund et 
al. 2006 Jl^). It has already been suggested (Cayrel et al. 2007 [|4|) that the convection-related 
asymmetry of the ^Li line could mimic the signal attributed so far to the weak blend of ^Li in the 
red wing of the ^Li line. However, this suggestion was only based on a hydrodynamical simu- 
lation for a single set of atmospheric parameters, representing the halo turn-off star HD 74000. 
Now the theoretical line asymmetry has been computed for an extended range in effective tem- 
perature, gravity and metallicity, covering the stars of the Asplund et al. sample. The computed 
asymmetry is about two per cent, a value which is also the mean of the ^Li/^Li ratio of the full 
sample determined with symmetric profiles by Asplund et al. 2006. Our conclusion is that these 
observations can be reinterpreted in terms of intrinsic line asymmetry, plus the amount of ^Li/^Li 
given by Asplund et al. reduced by 0.02. This drastically reduces the number of certain ^Li 
detections, from 9 to 1 or 2. 
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1. The WMAP / Spite plateau discrepancy 

The Big Bang nucleosynthesis leads now to very accurate predictions of the primordial abun- 
dances, with the cosmological parameter 77=A'^(baryons)/A'^(photons) constrained by the observa- 
tions of WMAP. On the usual scale of abundances by number of atoms per 10^^ hydrogen atoms, 
the logarithmic abundance of 'Li is log ?i(Li) = 2.64 ±0.04 (Spergel et al. 2007 [[l]]), against a 
logaritmic abundance of 2.22 it 0.08 for the Spite plateau (Charbonnel & Primas 2005 |^). The 
first idea which comes to mind is that the nuclear fragility of 'Li can deplete the element by the 
same mechanism which has depleted this nuclide in the Sun by about 2 dex in 4.6 Gyr. The main 
difficulty with this scheme is that the Spite plateau is both flat and shows only a very small scatter. 
It is a priori difficult to imagine a mechanism leading to a constant depletion over two decades in 
metallicity and almost 1000 K in effective temperature. The problem of the small scatter was also 
serious as long as effects of stellar rotation were believed to be a main ingredient of the depletion 
mechanism. The inital distribution of stellar rotation induces without doubt a statistical scatter. We 



shall come back to this mechanism in subsection 1.3 



Other possible mechanisms have been proposed. In section [O] we shortly discuss the proposal 
that the depletion of the Spite plateau is due to a massive astration of the cosmic matter before 
the birth of the galactic halo stars (Piau et al. 2006 [^). In section L2 we quote proposals using 



supersymmetric particles, very briefly, because the subject is considered elsewhere in this volume. 



In section 1.3 we discuss in detail the diffusion - slow mixing - nuclear burning mechanism, now 



widely studied in dwarfs of solar metallicity. 

1.1 A major astration of cosmic matter before tlie formation of halo stars 

Piau et al. 2006 [Q] have proposed that about half of the mass of the cosmic matter from which 
the halo stars formed was already processed by pregalactic supemovae (SNe) of masses between 
10 and 40 Mq. Although this scenario has some interesting features concerning, for example, the 
high depletion of Li in the extremely iron poor, C-enriched star HE 1327-2326, the level of astration 
invoked seems incompatible with the expected resulting metallicity on the part of the plateau below 
[Fe/H]=-2.0 (Prantzos 2006, NIC -IX g]). 

1.2 Supersymmetric particles 

If they exist, supersymmetric particles may have played a role, both in depleting 'Li and in 
producing ^Li during the Big Bang, winning on the two sides. As the topic has been developed 
elsewhere in this volume, I will just mention here a few references dealing with this topic: Jedamzik 
2004 Pospelov 2007 [^, Cumberbatch et al. 2007 [|], and Kusakabe et al. 2007 [|]. 

1.3 Depletion by the diffusion - slow mixing - nuclear burning mechanism 

Atmospheres of low mass unevolved stars are rapidly mixed by convection and a small amount 
of convective overshooting. Below the convection zone, the situation is rather complex. Heavy ions 
settle down very slowly, but if a new ionization stage is reached, radiative support by continuous 
plus line absorption can block this motion. Under the impulse of Georges Michaud (see the Meeting 



in Honor of Georges Michaud 2005 [10|) an enormous amount of work has been done in this field. 



The situation is complicated by the existence of hydrodynamical effects which tends to modify this 
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vertical stratification and by the existence of meridional circulation. Initially, turbulence generated 
by differential rotation was considered to be the main cause of mixing, inhibiting to some degree 
the effects of diffusive sedimentation. The presence of hydrodynamical turbulence which tends to 
scramble the effects of diffusion, also enables the depletion of fragile nuclides which are destroyed 
at higher temperatures some depth below the bottom of the convection zone. In that case, the 
hydrodynamical turbulence creates a slow mixing between the convection zone and the burning 
region below. Very clearly, this is the case for Sun, in which the photospheric ^Li is depleted by 
more than a factor 100 from an ISM initial value of 3.3. Korn et al. 2007 [|ll|] have constrained the 
role of diffusion by comparing the abundances of a few elements in unevolved stars and evolved 
stars in the globular cluster NGC 6397. They derive a depletion of ^Li by about 0.25 dex using 
a stellar model of O. Richard (see [10]). This value depends upon the difference of effective 
temperature between the unevolved stars and the evolved stars, and upon the empirical model of 
Richard, done before the theoretical work of Charbonnel & Talon 2008 [14]. The latter approach 
has been very successful in explaining a large body of observations able to constrain the theory. 
First the evolution with time can be checked by observing stars in clusters of various ages, from 
the Pleiades to M67. Then the Boesgaard-Tripicco dip [12] poses a severe challenge to the theory, 
which was resolved only very recently [14]. The simultaneous observation of the three fragile 
nuclides ^Li, ^Be, and ^B in field stars (Boeesgaard et al. 1998 [13]) supplies direct evidence that 
the depletion is linked to the physical properties of the nuclei, and reaches levels unattainable 
by diffusion only. This body of observations has allowed to put the theory of loss of angular 
momentum and of chemical mixing in cool stars on a much better footing than before. A key 
ingredient has been to take into consideration internal gravity waves as well as the effects of rotation 
(meridional circulation and shear turbulence) [14]. Unfortunately the low-metallicity stars of the 
galactic halo have not yet been investigated in as much detail as population I stars. Further work is 
needed before a well accepted conclusion can be drawn. 



1.4 New observational results on the metal-poor tail of the Spite plateau 

Very few unevolved stars at extremely low metallicity are bright enough to be subjected to the 
spectroscopic study of the weak Li feature. A sudden drop in the ^Li abundance at metallicities 
below 0.01 the solar metallicity is quite visible in a diagram plotted not in logarithmic coordinates, 
but on a linear scale (see Fig. 22 in Asplund et al. 2006 [|^], referred to as ALN06 further on). In a 



study of the unevolved metal-poor halo stars, Bonifacio et al. 2007 Q16| ] have studied 17 objects with 
metallicities between —3.8 < [Fe/H] < —2.7 (see Fig. [I]). This study supports the conclusion that 
in this metallicity range something is happening to the Spite plateau: (i) the scatter in Li abundance 
is much larger here than on the plateau (—2.5 < [Fe/H] < —1.5). The double-lined spectroscopic 
binary CS 22876-032 is a striking example of this behavior. Component A (higher mass) has a 
logarithmic ^Li abundance of 2.2, component B has a Li abundance of only 1.75 ([|T7|], see also 
Fig. [T]). (ii) a decline of the abundance of ^Li at very low metallicity looks probable. Garcia Perez 
2008 [ p^ also finds a drop at these lowest metallicities, which can be connected to the strict upper 
limit of ^Li < 0.6 found in the ultra-metal-poor star HE 1327-2326 at [Fe/H]= -5.3 (Frebel et al. 
2008 [III). 
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Figure 1: The metal-poor end of the Spite plateau. Note the increased scatter, and the trend towards a drop 
of the Li abundance at the lowest mean metallicities. Adapted from Bonifacio et al. 2007 [ p^ and Gonzalez 
Hernandez et al. 2008 [FtIi. 



2. The ^Li problem 

Until 2006, only very few determinations of the abundance of ^Li were performed. The ob- 
servation is very difficult, the resonance line of ^Li being a weak doublet, blended with the cor- 
responding doublet of ^Li, shifted by only 0.16 A to the red of the 'Li line, and about 20 times 
weaker. The doublets are not resolved because of the thermal and turbulent broadening of the 
hnes. So, the presence of ^Li is detected only by a slight extra depression of the red wing of the 
unresolved 'Li feature (Fig. ^. A drastic change occurred with the paper by Asplund et al. 2006 
[jlsl], who published a set of 24 determinations of ^Li abundances, based on high quality spectra 
acquired with ESO's VLT/UVES combination. The surprising result was a kind of ^Li plateau (top 
of Fig. 23 of their paper) at the level logM(^Li) w 0.8. As clearly demonstrated by their Figure 23 
this behavior does not fit the expectations based on cosmic ray (CR) production (again their Fig. 23, 
bottom). This result has immediately triggered new theoretical papers trying to explain the new ^Li 
plateau. 

2.1 The pregalactic Cosmic Ray solution 



Rollinde, Vangioni & Olive 2006 [ |20| ] have considered pregalactic cosmic rays generated by 
shocks around population III objects. By collisions between accelerated and at rest helium nuclei, 
the so called a fusion process generates both Lithium isotopes by the reactions '^He(a, D)^Li 
and '*He(oi, p)'Li. This process works, but needs a very large fraction of the kinetic energy of 



population III SNe to be converted into CR energy, as noted by Prantzos 2006 [ ]21[ ] and Evoli et 
al. 2008 [^]. It seems also difficult to produce a plateau - in the diagram log?i(^Li) vs. [Fe/H] - 
in the hierarchical model of Galaxy formation of Evoli et al. 2008. Other scenarios involving CR 
production have been proposed, but are exposed to similar critics. 

2.2 The supersymmetric particle solution 



See the list of references given in section 1.2, and paper by B. Fields (NIC-X, this volume). 
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Figure 2: left panel:Fig. 13 (middle) of ALN06. Selected objects: filled dots, discarded objects: open 
circles. The mean value of n(^Li)/«(^Li) of the full sample is 0.022, the mean value of the selected subsample 
is 0.041. Right panel: Same as left panel excepted tha the selection is now done on the accuracy only error 
< 0.013. The mean value is not affected now 



2.3 The ^Li asymmetry solution: new results 

Observationally, the ^Li measurement relies on the detection of a slight extra depression of 
the red wing of the 'Li feature. However, such a line asymmetry already exists without the pres- 



ence of ^Li. 3D hydrodynamical models of convection produced by the CO^BOLD code [23] 
clearly predict a convection-induced line asymmetry of similar strength. Cayrel et al. 2007 [24] 
have shown that, in the case of the star HD74000, a halo turn-off star of metallicity [Fe/H] = -2.0, 
the derived ^Li abundance is reduced to a negligible amount when taking into account this intrin- 
sic line asymmetry. It was also shown that iron lines with equivalent widths similar to those of 
the doublet components of the 'Li feature, and with the same stratification, show about the same 
amount of asymmetry as the computed one for 'Li. This is the only "sanity" check suggesting that 
our computed asymmetries are correct. The referee asked us if we could compare our asymmetries 
with those obtained with the 3D models used in ALN06, for example. This was not possible for 
'Li, by lack of published data. We tried to use the published data in Allende Prieto et al. [ ^5[ ] on 
iron bisectors for Procyon. However in the short time available we do not claim to have reached a 
definite conclusion. Work is in progress to investigate this important point. 

The big question is to know what amount of ^Li would remain in the Asplund et al. analysis 
of 24 halo stars if the natural asymmetry of 'Li had been properly included. In order to get some 
idea of that, we have extended the computation of the theoretical asymmetry of the 'Li feature to 
other values of the parameters metallicity, effective temperature and gravity, spanning the range 
of Asplund et al. observations. Specifically, new computations have been performed for the cen- 
tral effective temperature of the sample, reff = 6300 K, at the 3 metallicities [Fe/H]=— 3.0, —2.0, 
and —1.0, and for the two gravities logg = 4.0 (turn-off) and 4.5 (unevolved dwarfs). At the cen- 
tral metalliciy, [Fe/H]=— 2.0, the computations were made also for reff = 5900 K and again for 
logg = 4.0 and 4.5. The results of these computations are presented below. 
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2.3.1 A statistical problem with the subsample with detected ^Li 

Before proceeding further, it is worth to note that what everybody calls the ^Li plateau is 
the almost constant ^Li abundance of the stars having a certain detection of ^Li, namely those 
satisfying the condition: (^Li)/('Li) > 2a ; a being the estimated accuracy in the determination 
of this ratio. As the abundance of ^Li is almost constant in the sample, and the noise is also fairly 
constant, it means that all values of n{ ^Li)/n('Li) below about 2 per cent have been discarded. This 
is clear in Fig. 13 of ALN06, reproduced here as Fig. ^. We are interested in the meaning of the 
"plateau" with respect to the full sample we note that the mean of the ratio for the full sample is 
0.022, compared to 0.041 for the subsample. The median of the full sample is 0.018. A selection 
criterion involving the value of the studied variable is well known to introduce a severe bias (see the 
warning of Kapteyn & Weersma 1910 [26] about small parallaxes). A selection with ?i(^Li)/?i(^Li) 
> 3a would move the plateau to a still higher value of 0.045, whereas a selection on the error only, 
keeping errors below 0.013, when the mean error is 0.0154, leads to a mean of ?i(^Li)/?i(^Li)=0.021, 
almost unchanged (Fig. ^. A more severe constraint on the accuracy of the determinations, error 
< 0.012, gives a mean of 0.18, still close to the mean of the full sample. In what follows we claim 
that the true measured signal is about 2 per cent induced by the asymmetry of ^Li (affected by the 
measurement error). The mean of 0.022 acquires then a new signification: absence of ^Li. The 
certain detections are now decreased in number, as only values of 2a above 0.02 represent certain 
detections. The new number is of the order of 1 or 2, and critically depends on the value of a. 



2.3.2 The size and behavior of the convective line asymmetry 

Our main interest is investigate the dependence of the convective line asymmetry upon stellar 
metallicity, effective temperature and surface gravity. Table |I| gives our preliminary results for a 
small grid of 8 hydrodynamical simulations. The entry 'asym' is the excess of absorption in the red 
wing of the blend relative to the equivalent width of the full line. 'Asym' was determined by fitting 
the asymmetric 3D line profile, obtained from NLTE line formation calculations on a number (A^) 
of selected snapshots from the hydrodynamical simulations, with a combination of two symmetric 
ID line profiles. It is what is measured in practice to determine the ratio «(^Li)/w('Li). 

For given Teff and logg, the predicted line asymmetry is fairly constant for metallicities below 
[Fe/H]=— 2, but increases noticeably towards higher [Fe/H]. Hotter stars with lower gravity (turn- 
off stars) show higher line asymmetries than cooler unevolved dwarfs. 

2.4 The source of the line asymmetry 

It is interesting to have a look at the dynamical origin of the line asymmetry. Obviously, it 
is the hydrodynamical asymmetry between upward and downward convective flows that is at the 
origin of the asymmetry of the 3D line profiles. Fig. ^ shows a histogram of the vertical velocity at 
4 optical depths in a typical 3D hydrodynamical model (Tgff = 6250 K, logg = 4.0, [Fe/H]= —2.0). 
There is little asymmetry at log Tross = — 3 or — 2, but in the lower photosphere (log Tross = — 1 and 0) 
the asymmetry is conspicuous. For the centre of the disc, the 3D-NLTE contribution function of the 
Li I resonance line peaks at log Tross ~ — 1 . The broadening is large and the asymmetry is important 
(Fig. |5|). In 3D-LTE, however, the contribution function has a double peak, one at log Tross ~ — 1, 
but a second one at in very superficial and very cool layers [^], where the velocity distribution is 
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Figure 3: The red line is a (intrinsically asymmetric) 3D-NLTE profile of the major component of ^Li 
(area W), computed for atmospheric parameters Teff = 6250 K, \ogg = 4.0, [Fe/H] — —3.0. The blue line 
is the mirror image (with respect to the bisector of the line core) of the red line completing a symmetric 
profile with the left wing of the red profile. The difference between the right wings of the red and the 
blue profiles (area AW) measures the amount of asymmetry of the 3D profile (signal = AW/W — 0.020). 
The green dots represent the profile of a ^Li / ^Li blend, computed with the symmetric profile and a ratio 
«(^Li)/n(^Li)=0.022, equal to the mean value of the observations of ALN06 for a selection of stars with 
parameters 6100 K< Teff < 6400 K and 3.7 < \ogg < 4.3. The signal of this blend measures the isotopic 
ratio: «(''Li)/«(^Li) =AW/W =0.021. There is a slight difference between the signals produced by the 
asymmetric 3D profile and the ^Li / ^Li blend (the ^Li blend is slightly more red-shifted and shallower), but 
these differences are too small to be measurable. 
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Figure 4: Same as fig |[ except that here [Fe/H]= -2.0. 
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Figure 5: Histogram of vertical velocities at 4 optical depths in the atmosphere of a typical metal-poor star 
(Teff = 6250 K, \ogg — 4.0, [Fe/H]= —2.0). Positive (negative) velocities correspond to blue- (red-) shifts. 
See text for details 



Table 1: Computed line asymmetry expressed as equivalent isotopic ratio ^Li/^Li (asym.), and 3D model 
properties: Teff is the effective temperature, \ogg the surface gravity, [Fe/H] the metallicity, the number of 
considered statistically independent snapshots, and X xY xZ the geometrical extend of the computational 
box. All models use six wavelength bins to represent the wavelength dependence of the radiative energy 
exchange, and a grid of 140 x 140 x 150 cells. 

narrower and virtually symmetric. This explains why the 3D-LTE analysis of ALN06 differs from 
our 3D-NLTE analysis. There is a substantial difference in the two profiles, the 3D-LTE profile 
being a superposition of an asymmetric line (the part formed in the lower photosphere) and of a 
symmetric line (with a different shift) formed in the higher photospheric layers. 

Clearly, the real situation is more complex than suggested by the simplistic picture outlined 
above. It is not only the distribution of the vertical velocity but also its correlation with the temper- 
ature fluctuations that determines the shape of the line profile. Moreover, horizontal motions also 
become important when integrating the line profile across the stellar disk. 



Teff \ogg [Fe/H] asym. X x F x Z [Mm^] 



6250 4.0 -3.0 0.020 20 26 x 26 x 12.7 

6250 4.0 -2.0 0.021 16 26 x 26 x 12.7 

6250 4.0 -1.0 0.037 20 26x26x 12.7 

6250 4.5 -3.0 0.012 18 7.0x7.0x3.9 

6250 4.5 -2.0 0.011 19 7.0x7.0x3.9 

6250 4.5 -1.0 0.017 20 7.0x7.0x3.9 



5900 4.0 -2.0 0.015 20 26 x 26 x 12.5 
5900 4.5 -2.0 0.011 19 6.0x6.0x3.8 
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2.5 Conclusion 



The level of the plateau attributed to ^Li in ALN06 [15] is affected by a bias which, when cor- 
rected, moves the plateau to a mean value of ^Li/^Lif» 0.022. The study of the 'Li line asymmetry 
with a 3D-NLTE code on 3D hydrodynamical models obtained with the CO^BOLD code leads to 
asymmetry signals of respectively 0.020, 0.021 for [Fe/H]= —3.0 and —2.0 respectively, at Teff = 
6250 K, \ogg = 4.0. Stars corresponding to the subsample of ALN06 with 6100 K< Tgff < 6400 K, 
3.7 < logg < 4.3,have a ratio ?i(^Li)/?i('Li) with a mean value of 0.021 (14 objects). For metallicity 
[Fe/H]=— 1, the predicted asymmetry signal is 0.037, and the observed value in ALN06 is 0.023, 
not a very good fit but the number of objects is only 4. The effect of the gravity is a decrease of 
the asymmetry by a factor of 2 when logg = 4.5 instead of 4.0. The only certain unevolved star in 
the ALN06 sample is HD 19445: observed signal: 0.002, asymmetry predicted 0.010. The conclu- 
sion is that the plateau of the remarkable ALN06 set of observations can be reinterpreted without 
a contribution by ^Li. However, a detailed re-analysis of a few stars with the highest apparent ^Li 
content is necessary before a final conclusion can be drawn for them. 
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